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ABSTRACT 

Published mid infrared spectra of transiting planets HD 209458b and HD 189733b, obtained during sec- 
ondary eclipse by the InfraRed Spectrograph (IRS) aboard the Spitzer Space Telescope, are predominantly 
featureless. In particular these flux ratio spectra do not exhibit an expected feature arising from water vapor 
absorption short-ward of 10 /im. Here we suggest that, in the absence of flux variability, the spectral data for 
HD 189733b are inconsistent with 8 /im-photometry obtained with Spitzer' s InfraRed Array Camera (IRAC), 
perhaps an indication of problems with the challenging reduction of the IRS spectra. The IRAC point, along 
with previously published secondary eclipse photometry for HD 189733b, are in good agreement with a one- 
dimensional model of HD 189733b that clearly shows absorption due to water vapor in the emergent spectrum. 
We are not able to draw firm conclusions regarding the IRS data for HD 209458b, but spectra predicted by 
ID and 3D atmosphere models fit the data adequately, without adjustment of the water abundance or reliance 
on cloud opacity. We argue that the generally good agreement between model spectra and IRS spectra of 
brown dwarfs with atmospheric temperatures similar to these highly irradiated planets lends confidence in the 
modeling procedure. 

Subject headings: stars: planetary systems, individual (HD 189733, HD 209458), radiative transfer 



1. INTRODUCTION 

Since the detection of the first extr asolar giant planet (EGP), 
51 Peg b (Mavor & Oueloj Il995h . considerable effort has 
gone into both observing and modeling properties of the at- 
mospheres of the close-in "hot Jupiter" planets. The Spitzer 
Space Telescope now provides us a unique probe into the mid- 
infrared emission from these exotic atmospheres. 

In this Letter we examine mid infrared spectra of transit- 
ing planets HD 209458b and HD 189733b, which were re- 
ce ntly observed with Spitze r's In fraRed Spectro g raph ( IRS) 
by iRichardson et all (l2007l) and iGrillmair etall (l2007l) . re- 
spectively. The spectra were obtained from ~7.5 to 14 iim 
around the time of secondary eclipse. The most prominent 
conclusion of these works is that strong absorption arising 
from H2O, expected short-ward of 10 /im, was not seen for 
either planet. Possible explanations put forward for the rela- 
tively featureless spectra included a lack of H2O in these at- 
mospheres, masking by a continuum of opaque clouds, or an 
isothermal tem perature str ucture. 

Recently iKnutson et al.l (l2007al) published an 8 fim light 
curve for HD 189733b that includes an observation of the sec- 
ondary eclipse. These observations used Spitzer's InfraRed 
Array Camera (IRAC). We suggest that this data point is in- 
consistent with the IRS spectrum, that water vapor is present 
and absorbs from 6.5-10 /^m, and that the published, low S/N 
spectrum of HD 189733b may not reflect the true planet spec- 
trum. We have computed model atmospheres and spectra of 
HD 189733b and HD 209458b in an effort to shed additional 
light on these exciting observations. We will also briefly dis- 
cuss L- and T-type brown dwarfs, whose atmospheres span 
the same temperature range as hot Jupiters. We suggest that 
the generally very good fits of our models to these high signal- 
to-noise IRAC and IRS datasets supports the veracity of the 
underlying physics and chemistry of our planet models. 
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- Spitzer Fellow, Carl Sagan Center, SETI Institute, 515 North Whisman 
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2. SPITZER SPECTRA OF PLANETS 

Two teams employed the IRS instrument in "Short-Low" 
configuration to obtain low resolution spectra of parent stars 
HD 209458 and HD 189733 around the time of planetary sec- 
ondary eclipse. Each group compared the spectrum of the 
system when the planet's day side was still visible (before 
and after the secondary eclipse) to the stellar spectrum when 
the planet was behind the star. The observed enhancement 
in flux seen when the planet was visible was attributed to the 
planet, thus providing a measurement of the planet-to-star flux 
ratio. These exciting and difficult observations are at the lead- 
in g edge of spectr a l chara cterization of exoplanets. 

IGrillmair et al.l (120071) observed two secondary eclipse 
events of planet HD 189733b and from these observations 
published a planet-to-star flux ratio from 7.57 to 14.71 fim. 
They found a flux ratio spectrum that is essentially flat across 
these wavelengths with a mean eclipse depth of 0.0049 ± 
0.0002. The subset of their data with the highest signal-to- 
noise is shown in Figure[T}i. Richardson et al. (2007) ob- 
served two secondary eclipses of planet HD 209458b and pub- 
lished a flux ratio spectrum from 7.5 to 13.2 /im. The derived 
average ratio is shown in Figure [!}:. These authors also find a 
spectrum that could be characterized as mostly flat. However, 
they also detect increases in the flux ratio from 9.3-10 /im, as 
well as a narrow emission feature at 7.7 /im, which they at- 
tribute to emission from a stratosphere. 

In addition to spectra, we plot Spitzer secondary eclipse 
photometric data. These points include the IRS blue peak- 
up at 16 /im (Deming et al. 2006) and the IRAC 8 /im point 
(Knut son et al.ll2007al) for HD 189733b, and the Multiband 
Imaging Phot ometer for Spitzer ( MIPS) point at 24 /tm for 
HD 209458b dDeming et al.ll2005h . We note that in FigurefTJz 
the IRAC 8 pm point clearly plots below the IRS data, which 
we'll discuss in detail in §2.2. 

2.1. Model Atmospheres and Spectra 

In order to further the study of these observations and their 
constraints on the character of each planet's atmosphere we 
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Fig. 1 . — (a) Planet-to-star flux ratio for planet HD 189733b. Black diamonds are the IRS secondary eclipse data froml Grillmair et aTl 12007 !). The blue circles 
are secondary eclipse data from Knutson et al. 1 2007a) at 8 ^m and Deming et al. ( 2006) at 16 fj.m. All error bars are Icr. The dotted curves are transmission 
functions for these photometric bandpasses. In gray are flux ratios for blackbody planetary emission at 1400, 1200, and 1000 K. The red solid curve is a ID 
uniform day/night model. The red dash-dot curve is a model with a hotter day side. Green diamonds are computed ratios for the dash-dot model across the 8 and 
16 fim bands, (h) Shown are absorption cross-sections (left axis) for CO (green), H2O (blue), and CH4 (red). In black are the opacities (Q„fc,/particle radius, right 
axis) of 1 fj,m radius spheres of amorphous forsterite (dashed), amorphous enstatite (dash-dot), and c rystalline enstatite (solid), (c) Planet-to-star flux ratio for 
planet HD 209 45 8b. Black diamonds are the IRS secondary eclipse data from Richardson et^ t2007|). The red circle is the secondary eclipse observation from 
IDeming et al.l (2005) at 24 ^m. The dotted curve is the blue edge of the transmission function for this bandpass. The gray curves are the planetary blackbodies 
from above. The red curve is from the ID uniform day/night model. The blue curve is from the day side of the 3D model from, Fortney et al. ( 2006a). 

have computed model atmospheres and generated low reso- 
lution spectra. We employ a ID model atmosphere code that 
has been used for a variety of planetary and substellar ob jects. 
Recently it has been used for brown dwarfs jM arley et al. 
2001 'Saumon et al.' '20061) and EGPs dFortnev et al.. ' ,2005,. 
2006b; Marlev et al. 2007). It explicitly includes both inci- 
dent radiation from the parent star (if present) and thermal 
radiation from the planet's atmospher e. The radiative t rans- 
fer solving scheme was developed by iToon et aTl (1 19891) and 
was first applied to irradiated giant planet atmospheres by 
Marlev et al. (1999). We use the elemental abundance data of 
LoddersI (12003ft and compute chemical equilibrium composi- 
tions following Lodders & Fegley (2002, 2006) . Our opacity 
database is des cribed in [Freedman et al.l (l2007h . The stellar 
models are from 'Kurucz'('l993l) and the system parameters are 
froml\yinn etal . (2007) for HD 189733b and lKnutson et all 
(12007b') for HD 209458b. The planet models presented here 
are cloud-free. 

We plot model planet-to-star flux ratios in Figure [T] For 
HD 189733b we show two one-dimensional models. The 
solid red curve is for a day-side model that assumes all in- 
cident flux is efficiently redistributed over the entire planet, 
meaning that th e incident flux is multiplied by a factor of 1/4 



(I Appleby & Hog an 1984). Following co nventions adopted i n 
solar system atmosphere modeling (e.g. iMcKav et al.lll989h . 



here and in our previous papers (iFortnev et al.ll200l l2006bl) 
the incident flux is decreased by incorporating a multiplicative 
factor / of 1/2 due to the day-night average, and another fac- 
tor of 1/2 arises from the mean 60° stellar zenith angle (/i, the 
cosine of this angle, i s 0.5) that is correct for global average 
insolation conditions (iMcKay et al.lll989t iMarley & McKayl 
Il999l) . With this prescription more stellar flux is absorbed 
at low pressures, leading to a warmer upper atmosphere and 
shallower temperature gradient, compared to a model where 
incident flux is directly multiplied by 1/4 with ^=1.0. The 
greater absorption of incident flux at lower pressure is due 
to the twice longer path length for incident photons to reach 
a given pressure. In reality, these differences in ID models 
could be swamped by a more complex day-side temperature 
structure for these planets. Here, for models with a warmer 
day side, the incident flux is multiplied by a factor / be- 
tween 1/2 and 1, while /i remains 0.5. The dash-dot curve 
in Figure[T}: is for a warmer model, where / = 4/5, meaning 
that only 20% of the absorbed incident flux is assumed to be 
advected to the night side. Both models predict an absorption 
feature due to H2O short- ward of 10 fim. Figure [27 shows ab- 
sorption cross sections for H2O, CO, and CH4. At these tem- 
peratures (^1200 K), CO is strongly favored at the expense 
of CH4. The dashed model in Figure [Th is a goo d (~lcr) fit 
to the band-averaged 8 fim dKnutson et al.ll20b7ah and 16 fim 
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dPeming et a l.|2006| ) secondary ecli pse observations. We note 
that although Grillm air et alJ (l2007h claim their flat ratio spec- 
trum is consistent with a blackbody, Figure [Hz shows that it is 
actually bluer than a blackbod^; 

Figure[TJ: shows the mean lRichardson et alJ (1200 7*) data set. 
In solid red we plo t the / = 1/2 ID model spectrum from 
iFortnev et alj (l2005l) . Note that the predicted "downturn" in 
the ratio spectrum short- ward of 10 /im is actually quite small 
in this ID mo del, less than that predicted by o ther model- 
ers (e.g.. lSeager et al .' 2005; Burro ws et aLll20 06'), for reasons 
that are not yet clear, but could be related to the angle of inci- 
dent radiation. Given the S/N of this dataset and the relatively 
small predicted depth of the model H2O absorption feature, it 
is not clear that absorption due to H2O can be dismissed. 

Dynamical processes can act to produce isothermal profiles 
in the upper atmosphere which would also suppress the depths 
of absorption feat ures. In sohd blue we plot t he spectrum of 
the day side of the lCooper & ShowmanI (l2006h 3D dynamica l 
model of HD 209458b. as computed in lFortney et a!I(l2006ah . 
The isothermal atmosphere is due to a super-rotating west-to- 
east equatorial jet; as cool night side gas is blown around the 
planet's limb this gas is more readily warmed at lo wer pres- 
sures where the radiative time constant is shorter (liro et alJ 
l2005h . leading to day side P-T profiles that are more isother- 
mal than one would predict for a radiative equilibrium model 
ISh owman & Cooper 2006). 

Neither the ID nor 3D model has the increase in flux ra- 
tio seen from ^9.3- 10 /im, as neither model has the strong 
temperature inversion that would be necessary to create emis- 
sion features. In Figure[Tf' we plot the Mie absorption effi- 
ciency of 1 /im silicate particles of forsterite (Mg2Si04) and 
enstatite (MgSi03). Larger particles have a gray er spectrum 
and d o not exhibit the canonical silicate feature (.Hanner et alJ 
II994I) . The location of this possible emission feature in the 
HD 209458b spectrum is suggestive, but it will take higher 
signal-to-noise observations to determine if silicates have in- 
deed been detected. HD 209458b is a better candidate than 
HD 189733b for having silicate clouds in its visible atmo- 
sphere, due to the higher stellar insolation and hotter atmo- 
sphere. 

2.2. A Closer Look at HD 189733b 

As was noted abo ve and shown in Figure fT};, the IRAC ob- 
servation at 8 /xm by lKnutson et al.l (l2007al) indicates a down- 
turn that is in fact consistent with H2O vapor absorption. We 
now take a closer look in Figure |2] Our model (re^)js a good 
fit to the data, while i n blue we show alBurrows et all ( l2006i) 
model, as taken from iGrillmair et aP (120071) . that appears to 
also be a good fit to the photometry, albeit with a much deeper 
water feature. Binning the 31 IRS data points from 7.57 to 
9.39 /im (which overlaps the IRAC band) into a single point, 
we find an eclipse depth of 0.00430 ± 0.00017, -5.3cr re- 
moved from the Knutson et al. IRAC point. If variation in 
plan etary or stellar flux were to account for this difference 
(e.g.' Menou et alJl2003HRauscher et al.ll2007!) . the variabiUty 
would have to be ~25% near 8 /im. If there is no variability 
and both the IRS spectrum and IRAC data point are accurate 
as published, we can examine how small the flux ratio must 
be short-ward of the bluest IRS point at 7.57 /xm. We have 
created synthetic data with the same spectral resolution of the 
IRS data set. In triangles we show a constant flux ratio of 
0.00156 between 6.3 //m (the blue end of the IRAC filter) and 
7.57 /im. In boxes we show a synthetic ratio that decreases 
linearly from the 0.00430 binned value. The ratio must fall 
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Fig. 2. — Planet- to-stai' flux ratio for planet HD 189733b. In diamonds 
are the lGrillmair et a l. (2007) IRS data. The blue circles are from secondary 
eclipse photometry. The scaled band transmission functions are shown as 
dotted lines. The black circle is the eclipse depth from 7.57 to 9.39 ^m from 
the IRS data. All error bars are Icr. The triangles and squares show flux ratios 
short-ward of the IRS data that are small enough to bring the flux ratio across 
the entire IRAC band (6.3 to 9.4 /tm) into precise agreement with the IRAC 
observat ion. (See te xt.) In solid red is our / = 4/5 model and dash-dot blue 
is the Bu rrows et aU (2006) model as shown in Grillmair et al. (2007). 

to zero by 6.79 /im to give a band-average eclipse depth that 
agrees precisely with the IRAC data point. 

In Figure (TJ? we showed that the only potential absorbers 
predicted from equilibrium chemistry from 6 to 10 /xm are 
H2O and CH4, neither of which possess a sharp increase in 
opacity short-ward of 7.7 /im. With the assumption that the 
IRAC secondary eclipse depth is accurate, and given the se- 
vere absorption feature that would be needed just short of 
7.57 jjm to fit this point given the IRS spectrum, we conclude 
that the published IRS spectrum short-ward of 10 /im does 
not reflect the true spectrum at these wavelengths, where the 
IRAC data point suggests absorption by water vapor 

One possibility toward reconciling the data sets would be 
to shift the e ntire IRS spectrum d ownward. Indeed, for 
HD 209458b, iRichardson etaTI (120071) were unable to accu- 
rately determine the absolute depth of the secondary eclipse 
with the IRS, and therefore tied their spectrum to an un- 
published IRAC 8 /im observation of the sec ondary eclipse. 
Given this issue, the absolute depths of the IGrillmair et al.l 
(120071) data may also require a shift. However, the down- 
ward shift that would be required (^^25%) woul d lead to much 
poorer agreement between this data set and the lDeming et al.l 
d2006) point at 16 /im. In addition, although both teams 
corrected for slit losses, this was only done in a gray man- 
ner Since the point spread function is larger at longer wave- 
lengths, it is possible that slit losses would be more pro- 
nounced at the red end of the spectrum. Greater flux at longer 
wavelengths, together with a downward scaling, may be able 
to bring the observed spectrum for HD 189733b into agree- 
ment with the photometric points. Future work on all of 
these data sets may be beneficial, as, for example, the rela- 
tive photometry for these measurements is more than an order 
of magnitude mor e precise than IRAC's initial specifications 
(iFazio et al.ll2004 . 

3. LESSONS FROM BROWN DWARFS 

The P-T conditions encountered by the atmospheres of hot 
Jupiters is coincident with that found in the atmospheres of 
L- and T-type brown dwarfs. Although the actual P-T profiles 
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of irradiated EGPs are of course different, the underlying at- 
mospheric physics and chemistry are similar and indeed the 
actual computational models are in some cases virtually iden- 
tical, such that a comparison of brown dwarf models to data 
can illuminate theoretical areas that are well understood. 

Comparisons of our group's models to near- and mid- 
infrared spectra of L- and T-dwarfs, including superb observa- 
tions with the IRS, have generally found very good to excel- 
lent matches with data (Golimowski et al. 2004; Knapp et al. 



2004; Gushing et al. 2006; iSaumon et alj|200 6; Mainze r et alj 
2007; Leg gett et a l. 2007). S pecifically our "clou dy" L-dwarf 
and "clear" T dwarf models (iMarley et al.ll2002l M. S. Mar- 
ley et al. in prep.) fit the depths of the water bands lying 
betwe en Y, J, H, and K bands as well as the 6.5 Atm water 
band ( Gushing et al. 2006; Saumon et al."2006; Mainz er et alj 
[2007; Gushing et al. 2007) expected in the IRS spectra of hot 
Jupiters. There are, however, departures from the data in re- 
gions where molecular opacities are poorly known as well as 
some systematic spectral shape differences in K and N bands. 

Since the Marley et al. models have been demonstrated to fit 
the near and mid infrared data of the clou d-free T dwarfs w e 
gain confidence i n our radiative transfe r (iToon et al.lll989l) . 
opacity database ('Freedman et al."2007'), and chemistry cal- 
culations dLodd ers & Feglev 2002, 2006). This experience 
gives us confidence these areas are becoming well understood. 
Of course, since isolated brown dwarfs lack incident flux that 
drive atmospheric dynamics and photochemistry, these addi- 
tional areas may remain challenging topics of hot Jupiter char- 
acterization for some time. It is in deviations from ID or sim- 
ple day/night models that these aspects will become apparent. 

4. DISCUSSION & CONCLUSIONS 

The common conclus ion of the Grillmarr et al.l (l2007h and 
iRichardson et al.l (l2007h secondary eclipse observations was 



the apparent lack of absorption due to water vapor However, 
for HD 189733b we have shown that photometry indicates a 
planet-to-star flux ratio that is not as featureless as indicated 
by the published IRS data. If variability in planetary flux is to 
make up for this discrepancy, it must reach ^-^25%. The appar- 
ent downturn in the flux ratio from 16 to 8 fim implied by pho- 
tometry is consisent with a model that predicts absorption due 
to water shortward of 10 /im. As was shown in Figure |2] there 
is considerable diversity in feature depths among models, per- 
haps due to differences in chemistry, opacities, or treatment of 
incoming stellar flux. 

For HD 189733b, here we have shown that our model is 
a good match to 8 and 16 /im-band photometry. In addition, 
both our ID and 3D models fit the HD 209458b observations 
reasonably well. Our conclusion that water vapor absorption 
is present in th ese atmosphere s is sup ported by the recent 
analysis of the iKnutson et al.l (l2007bl) transit data set by 
iBarmanI ( l2007h . who concluded that water absorption in 
present in the atmosphere of HD 209458b in a band at ~ 950 
nm. Of course the characterization of the atmospheres of 
these planets is just beginning. Both planets have been 
observed in all 4 IRAG bands, but this data has not yet 
been published. Additional IRS spectroscopic observa- 
tions for HD 209458b (Deming and collaborators) and 
HD 189733b (Grillmair and collaborators) will be obtained 
this year. Importantly, these observations will go down to 
5.5 /im, fully overlapping the IRAG 8 fim band, which will 
allow for more accurate normalization of the obtained spectra. 

We thank J. Richardson, M. Gushing, D. Deming, G. Grill- 
mair, D. Gharbonneau, and the referee for useful discussions. 
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